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We are interested in finding the joint distribution function of the real and imaginary parts of the local Green 
function for a system with chaotic internal wave scattering and a uniform energy loss (absorption). For a 
microwave cavity attached to a single-mode antenna the same quantity has a meaning of the complex cavity 
impedance. Using the random matrix approach, we relate its statistics to that of the reflection coefficient and 
scattering phase and provide exact distributions for systems with (3=2 and /3=4 symmetry class. In the case 
of P=l we provide an interpolation formula which incorporates all known limiting cases and fits excellently 
available experimental data as well as diverse numeric tests. 
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Characterising statistical fluctuations of physical observ- 
ables in quantum systems with underlying chaotic classical 
dynamics remains a very active field of research in theoretical 
and experimental physics. A considerable progress in under- 
standing the phenomenon was underpinned by revealing the 
apparent universality of the fluctuations in systems of very di- 
verse microscopic nature, ranging from atomic nuclei and Ry- 
dberg atoms in strong external fields, to complex molecules, 
quantum dots, and mesoscopic samples, see e.g. |JJ|. From 
theoretical side, the universality allows one to exploit the ran- 
dom matrix theory (RMT) as a powerful tool for analysis of 
generic features of the energy spectra of such systems |2|,|3|]. 

In many atomic, molecular, and mesoscopic systems the 
quantity which is readily obtained experimentally is the ab- 
sorption spectrum for transitions from a given initial state 
\g) to highly excited chaotic states at the energy E. For 
high-resolution experiments chaotic spectra consists of well- 
resolved narrow resonance peaks and one can, in principle, 
study statistics of the peak heights and widths as well as that of 
spacings between consecutive peaks. Most frequently, how- 
ever, the absorption spectra look practically continuous due to 
both inevitable level broadening and finite experimental res- 
olution. Then the relevant statistics is the distribution and 
correlation functions of the absorption probability a(E) (also 
known as the strength function of the dipole operator fi) which 
in the simplest situation of uniform level broadening T is 
given by a(E) oc lm(g\fi G(E) fi\g), see e.g. Ref. |4] and dis- 
cussion therein. Assuming the validity of the RMT, the prob- 
lem then amounts to studying statistical properties of the re- 
solvent (Green function) operator G(E) = (E + iT/2-Hy 1 
associated with the random matrix H which replaces the ac- 
tual chaotic Hamiltonian. In particular, the imaginary part of 
the diagonal entries of G(E) is well known in the solid state 
physics as the local density of states (LDoS) and in this ca- 
pacity its statistics enjoyed many studies ||4j,|5l|6lLZllSll- 

From the experimental point of view the same universality, 
which makes the use of the RMT legitimate, provides one with 
an attractive possibility to employ simple model systems for 
analyzing generic statistics of the fluctuating quantities. One 



of such systems which proved to be an ideal playground for 
investigating a variety of quantum chaos phenomena are vari- 
ous microwave billiards 1 1]. The billiards are realized as res- 
onators in a form of electromagnetic cavities fed up with the 
radiation of the wavelength small compared with the charac- 
teristic size of the enclosure and shaped to ensure the chaotic- 
ity of internal scattering. The cavities are coupled to trans- 
mission lines or to waveguides which are used to inject elec- 
tromagnetic waves into the system as well as to collect the 
outgoing waves. The adequate description of the experiment 
is then achieved in terms of the scattering matrix S relating 
amplitudes of incoming and outgoing waves. 

High-resolution experiments are usually performed in low- 
temperature (superconducting) cavities with very high quality 
factor 12|]. Majority of the experiments are, however, done at 
room temperatures IKJ [TIJ [ljl [lj] . Inevitable energy losses 
(absorption) leading to uniform broadening of the resonances 
play, therefore, an important role and have to be taken into 
account adequately when describing the experiments theoret- 
ically. 

In the present paper we are going to concentrate on the sim- 
plest case of a single one-channel antenna experiment. To be 
able to employ the RMT methods, it is conventional to repre- 
sent the resonant part of the scattering matrix S in the follow- 
ing form (see e.g. Ref. J3): 



S(E) 



1 - iK(E) 
l + iK(E) 



ire 



(1) 



where K(E) = V^G(E)V is the so-called K matrix. The 
Hamiltonian H of the closed chaotic cavity gives rise to N 
eigenfrequencies E n characterized in the relevant range of 
the scattering energy E by the mean level spacing A. The 
column-vector V describes the energy-independent ampli- 
tudes coupling the corresponding eigenmodes to the propagat- 
ing mode in the antenna. We see again that the study of sta- 
tistical properties of the scattering matrix amounts to knowing 
statistics of the diagonal element G\i of the Green function of 
the closed system in some basis. In fact, in the present con- 
text the function iK = Z has the direct physical meaning of 
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the electric impedance Z of the cavity which relates linearly 
voltages and currents at the antenna port II ill . 

Without absorption T = and the scattering matrix is uni- 
modular: r = 1. At finite absorption, T > 0, the reflec- 
tion coefficient r and the scattering phase 9 have nontrivial 
distributions, which have been recently measured in experi- 
ment l id 1 1 3fl . On the other hand, universal fluctuations of 
both real and imaginary parts of the cavity impedance Z have 
been recently investigated experimentally in 1 11]. Since the 
impedance matrix Z is related to eigenmodes and eigenfre- 
quencies of the closed cavity, the study of Z is in some sense 
complementary to that of S. 

Various statistics related to the scattering matrix of chaotic 
systems with losses were subject of a number of recent papers 
O HIS El El 0. Explicit analytical results were avail- 
able, however, only for the simplest case of systems with no 
time-reversal invariance (TRI) corresponding to the so-called 
(3=2 symmetry class of the RMT. At the same time, major- 
ity of the billiard-type experiments is performed in systems 
which are time-reversal invariant ((3=1 symmetry class of the 
RMT). Similar situation holds for the statistics of the local 
Green function, in particular, LDoS. For the 0=2 case the cor- 
responding expressions were obtained by various methods in 
01 00 EH An attempt |8| to provide an expression for the 
LDoS distribution for the (3=1 case can not be considered as 
particularly successful, as the general expression was given in 
a form of an intractable fivefold integral. Finally, it is worth 
mentioning the existence of the (3=4 symmetry class describ- 
ing time-reversal invariant chaotic systems with half-integer 
spin. This situation may occur in quantum dots with strong 
spin-orbit scattering 1 22], in Rydberg atoms driven by mi- 
crowave fields 1 23], and can be efficiently simulated in some 
other models of quantum chaotic systems 13 • 

The fundamental quantity which determines the full statis- 
tics of S or Z is the joint distribution function V(u, v) of the 
real u = ReK and imaginary v = — ImK > parts of K. 
Generally, one can always write K = k(NA/tt)Gh in the 
RMT. The effective coupling constant k = w\\V\\ 2 /NA > 
enters the S matrix statistics only through the so-called trans- 
mission coefficient T = 1 — |S| 2 (= 4k / ( 1 + k) 2 in the middle 
of the spectrum, E = 0), see e.g. Ref. 1 14] for details. 

A convenient starting point of our analysis is the observa- 
tion that the distribution V(u, v) must always have the follow- 
ing general form: 



V{u,v) 



1 



2ttv' 



(2) 



with x = (u 2 + v 2 + l)/2v > 1, It initially emerged in J2J 
in the course of explicit calculations for the (3=2 symmetry 
class, but neither origin nor generality of such a form were 
appreciated. Here we show that Eq. (|2j is the direct conse- 
quence of two fundamental properties of the so-called "per- 
fect coupling" case T=l: (i) the statistical independence of 
the S matrix modulus r = |+j and its phase 9; and (ii) the 
uniform distribution of 9 e (0, 2ir). Both these properties can 
be verified using the methods of Ref. 112611 . The joint distri- 



bution P(x, 9) then factorizes to Pq(x)/2it. Choose now new 
variable y = ^| = cot(6l) instead of 9, so that \d9/dy\ = 
(l+y 2 ) -1 . Noticing that y = (u 2 +v 2 — l)/2u and evaluating 
the corresponding Jacobian \d(x, y)/d(u, v)\ = (1 + y 2 )/v 2 
of the transformation from x and y to u and v, we come after 
a simple calculation to (0. 

The explicit form of Pq(x) at arbitrary absorption for vari- 
ous symmetry classes will be given and discussed below. Hav- 
ing Pq (x) at our disposal, it is immediate to find the distribu- 
tion of the imaginary part v (the LDoS normalized for conve- 
nience to have the unit mean value): 



V v (v) = 



V2 

TV 3 / 2 







/ dqP Q 




h 





(3) 



The distribution is normalized to 1 and has the first moment 
unity, (v) = L dvvV{v) = 1, automatically due to invari- 
ance of the integrand with respect to the change v — > 1/v. 
Similarly, one can find the distribution of the real part u to be: 



V u {u) 



1 



27rVu 2 + 1 Jo 



dqP Q 



x/u 2 + 1 / 1 
—2—{ q+ q 



(4) 



Although u has no direct physical meaning in the context of 
solid state mesoscopic systems, both V u (u) and V v (v) are di- 
rectly measurable in microwave cavities II ill . 

Let us now discuss the explicit forms of Po(x) for various 
symmetry classes: (3=1, 2, 4. For the simplest case of broken 
TRI ((3=2) an exact result is available at arbitrary values of the 
dimensionless absorption strength 7 ee 27rr/A (3 03 En- 
sealing the absorption parameter for the subsequent use as 
a ee 7/3/2, we can represent the (3=2 result as follows: 

«M - f 

with a-dependent constants A = e a — 1 and B = 1 + a — e a 
and the normalization constant Mi = 1. 

For the case /3=4 the exact form became available very re- 
cently 1 3 1 ] by exploiting important advances in the RMtEI. 
The explicit derivation will be given elsewhere, the final result 
being 



A (a(x + l)/2) 



P/2 



D 



-a(x+l)/2 



(5) 



<piici 
131]: 



P gu »+C(x,7)e- 



7 sinh t 
dt 



(6) 



where P^ e (x) is the distribution (0 for (3=2 taken, however, 
at a = 27 and C(x, 7) ee j 2 (x+ l) 2 /2 - 7(7+ 1)0+ 1) +7. 

Unfortunately, for the most interesting case /3=1 the ex- 
plicit formula for Po(x) is not available yet, apart from the 
limiting cases of weak or strong absorption: 



Po(x) 



2r(/3/2+l) 
ae -a(«-l)/2 



(x+l)/2 



The first line here results from the relation 



7 < 1 

7 > 1. 
= 1- r 



(7) 



rl - - - ~ IT 

at 7 1 between the reflection coefficient and the (dimen- 
sionless) time-delay t, see flU Wfa . The time-delay dis- 
tributions are known Q Elfbr all (3=1,2, 4: V t (t) = 
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FIG. 1: The distribution Pq(x), Eq. ^5), at different values of the 
absorption parameter a=/3nT / A. Solid lines correspond to the 
(3=1 (GOE) and 2(GUE) cases and to the exact /3=4(GSE) result, 
Eq. Symbols stand for the numerics done for 400 realizations of 
500x500 random matrices. 



semble. The overall agreement of Eq. @ with the data at 
(3=1 is nearly as good as for the exact cases (3=2, 4. Another 
check has been performed in Ref. 1 13] which measured dis- 
tributions of the reflection coefficient and the scattering phase 
in the broad range of system parameters and found very good 
agreement with the corresponding expressions following from 
Eq. 0. 

We discuss now the behavior of V v (v) and "P u (u). Per- 
forming the integration in Eq. Q, we arrive for (3=2 at 

V°f e (v) = ( 7 /16^) 1 / 2 i;- 3 / 2 exp[-7(i; + i;- 1 )/4] 



x [2 cosh - + (v + v 



2/ 7 )sinh^], (8) 

which is exactly the LDoS distribution obtained earlier in @|. 
For the case (3=1, we get the following result: 



7>*>) 



27 v 3 / 2 



(A[K (a) + K t (a)]a + <firBe 



(9) 

with a = jq{\/v + 1/s/v) 2 and K v {z) being the MacDonald 
function. It is instructive to consider the asymptotic behavior 
of these functions in the limits of small or large absorption. 
At 7 <C 1, the distribution V v (v) becomes very broad having 
the maximum at v ~ 7, a power-law bulk behavior and an 
exponential cutoff at the far tails: 

, a (l+f3)/2 v -(3+P)/2 e - a /4v ^ y <<; a 

V v (v) cx < a 1 / 2 w- 3 / 2 , a< u< 1/a (10) 

{ a (l+/3)/2 w -(3-«/2 e -a,/4 > 1 / a ^ V] 

where constants ^ 1 are omitted. This result can be physically 
interpreted in the single-level approximation |5|,|8j], when the 
bulk and tail behavior is governed by spectral and wave func- 
tion fluctuations, respectively. As 7 increases, a number of 
levels contributing to v grows as ~ 7, so that V v (v) tends to 
the limiting Gaussian distribution 



V v (v) 



4 7ru 3 



exp 



a 

"I'v' 



1 ^~ 
7* 



(11) 



[(/3/2)' 3 / 2 /r(/3/2)]T- / 2 - 2 e- /3 / 2T . In the opposite case 7 > 
1, the known limiting Rayleigh distribution 11611 P{r) ~ 
(7/3/2)e- 1 " 1 '' 9 / 2 yields the second line in Eq.0 

In the absence of the general expression for (3=1 a natu- 
ral idea is to try to invent a formula interpolating between the 
known limiting cases 1 13]. We suggest Eq. to be the ap- 
propriate natural candidate, with the normalization constant 
being Afp = a/(AT(/3/2+l, a) + Be~ a ), where T(v, a) = 
f^°dtt IJ ~ 1 e~ t . Indeed, such a form reproduces correctly 
Eq. as both limits are determined solely by the first (uni- 
versal) term in (|5J- One needs, however, to keep B in order to 
handle properly the case of moderate absorption (a ~ 1). 

Figure 1 shows results of numerical simulations with ran- 
dom matrices drawn from the Gaussian orthogonal (GOE, 
(3=1), unitary (GUE, (3=2) or symplectic (GSE, (3=4) en- 



which has a peak at v ~ 1 of the width cx 1/^Jj 1, in 
agreement with the earlier result |8|] . 

As to the distribution V u (u), equation @ leads after the 
integration to the following exact (3 = 2 result: 



vriu) = j- 

Z7T 



2 °\ 2 



coshf /twn 
u X \2 ) 



(12) 

where u = y/u 2 + 1. Integrating the interpolation formula for 
the case (3 = 1, we obtain 



K oe {u) 



2-ku 



A py 



D 



(!) 



BKi ( 



4 J 



(13) 

where D(z)= f™dqy/l + z{q + q-^e'^+i ^ 4 is intro- 
duced for convenience. The limiting forms of V u (u) at weak 
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and strong absorption follow readily as 

7T _1 (1 + m 2 ) -1 , a<l 



V u {u) 



-au 2 /i 



(14) 



, a > 1 , 



and describe a crossover from the Lorentzian to Gaussian dis- 
tribution as absorption grows. This type of behavior as well as 
the trend of V v {v) to the Gaussian dl It was recently observed 
in the experimental study of the cavity impedance 111 ill . 

Finally, let us mention that the case of nonperfect coupling, 
T < 1, can be mapped |28, 29] onto that of perfect one mak- 
ing use of the parametrization 1 30 ] So = (S — — T)/(l — 
— TS). Here So is the scattering matrix of the system in 
the perfect coupling case. Now x and 9 do correlate and, af- 
ter the evaluation in parametrization Q of the corresponding 
Jacobian, the joint distribution P(x, 9) is found to be 



P{x, 9) = ±-P (xg - V(z 2 -l)( 5 2 -i; 



COS ( 



(15) 



with g = 2/T — 1. Complementary to Eq. (0, equation 
dl5> provides an access to scattering observables. The inte- 
gration there over x yields the scattering phase distribution. 
In particular, when absorption vanishes, x — > oo (r — ► 1) and 
P (x) -> 5(l/x), giving readily P{x,9) = p{9)5{l/x), with 
the phase density p(9) = [2ir(g — yj g 2 — 1 cos 9)]~ 1 found 
earlier 1 29] . As another example, the distribution of the reflec- 
tion coefficient in terms of Pq (x) is given at arbitrary coupling 
by (cf. Eq.(5) in 11911 and see also I13lll8ll in this respect): 



Pr(r) 



2tt 



d9ir- 



(1-r)' 



2(g-^/g T =l,/¥cos6) 



1 



9 



(16) 

In conclusion, although rigorous analytical treatment of the 
(3=1 case remains a theoretical challenge it is worth stressing, 
however, that the suggested interpolation formulas should be 
sufficient for the most of practical purposes of comparison to 
the experimental/numerical data. Moreover, all physically in- 
teresting limiting cases, as e.g. Eqs. dlOi . dl 11 and dl4> . are 
already reproduced from the exact limiting statistics 0. At 
the same time, it is important to understand that an excellent 
performance of the interpolation formula (|5jl for (3=1 is noth- 
ing else as a lucky coincidence. Indeed, applying the same 
formula for the (3=4 case, we have found that apart from well 
reproduced limits of weak and strong absorption, an agree- 
ment with numerics at intermediate values j ~ 1 turns out to 
be by far not as good as in the GOE case. 
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